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TECHNICAL NOTE 5072 


A THEORETICAL INVESTIGATION OF THE AERODYNAMICS OF 

WING-TAEL COMBINATIONS PERFORMING TIME -DEPENDENT 

MOTIONS AT SUPERSONIC SPEEDS 

By John C. Martin, Margaret S. Diederich, 
and Percy J. Bobbitt 


SUMMARY 


A theoretical investigation is presented of the contribution of 
horizontal tails to the lift and pitching moment due to angle of attack, 
a constant rate of pitch, and a constant vertical acceleration. Numer- 
ical values of the aerodynamic coefficients associated with these motions- 
fire presented for a number of two-dimensional wing-tail combinations, a 
triangular wing-tail combination, and a number of rectangular -wing- 
triangular -tail combinations. 

Methods for calculating the flow fields behind wings with constant 
vertical acceleration are developed. Calculated results are presented 
for the upwash behind two-dimensional wings and for certain regions 
behind triangular and rectangular wings for a constant rate of pitch and 
for constant vertical accelerations. A method of treating unsteady aero- 
dynamics based on an infinite series of stability derivatives of succes- 
sively higher order is also presented. 


INTRODUCTION , 


The development of the linearized theory of supersonic flow has 
permitted a first-order evaluation of a number of stability derivatives 
for a wide variety of isolated wings. (For example, see refs. 1 to 9 *) 
The linearized theory may also be used in the calculation of 'the aero- 
dynamic derivatives of the wing-tail combinations. Considerations of 
the effects of the horizontal tails located behind wings of necessity 
entail a knowledge of the upwash induced by the wing. The upwash 
from wings performing steady motions can be calculated by use of the 
various methods presented in references 10 to 15 . A method of calcu- 
lating the upwash in the region of the viscous wake behind wings with 
a constant vertical acceleration at supersonic speeds is presented in 
reference 1 6 . Except for the methods and results presented in refer- 
ences l6 to 19 , little consideration has been given to the theoretical 
calculation of the contribution of the horizontal tail to the aerodynamic 
derivatives. 
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The primary object of the present paper is a theoretical investi- 
gation of the contribution of the horizontal tail to the pitching moment 
due to a constant vertical acceleration. An investigation of this 
pitching moment entails a knowledge of a number of other factors which 
in th ems elves are useful in other supersonic flow problems. Inasmuch as 
a knowledge of these other factors is required in the primary investi- 
gation, the present paper has a number of secondary objectives such as: 
(l) the establishment of a method of treating unsteady aerodynamics of 
aircraft by the use of an infinite series of stability derivatives of 
successively higher order, (2) the development of theoretical methods 
for the calculations of sidewash and upwash behind wings which have 
local angle -of -attack distributions which vary linearly with time, 

( 3 ) an investigation of exact and approximate methods for the calcula- 
tion of sidewash and upwash behind wings with constant vertical accel- 
eration, and (L) an investigation of exact and approximate methods for 
the calculations of the lift and pitching moment due to a constant angle 
of attack and a constant rate of pitch. 

The upwash behind two-dimensional wings with a constant vertical 
acceleration is determined and calculated results are presented. The 
upwash along the center line of the wake behind triangular wings with 
subsonic leading edges is determined for a constant rate of pitch and 
for a constant vertical acceleration; calculated results are presented 
for a number of triangular wings and Mach numbers. The upwash in the 
wake in the plane of the wing behind rectangular wings is determined for 
a constant vertical acceleration; calculated results are presented for a 
number of rectangular wings and Mach numbers. 

Exact linearized results are presented for the lift and pitching 
moment resulting from angle of attack, steady pitching, and a constant 
vertical acceleration for a number of two-dimensional wing-tail combina- 
tions, a triangular wing -tail combination, and a number of rectangular- 
wing — triangular -tail combinations. Some of these exact results are 
compared with results calculated from a number of simple approximate 
relations which are often used in calculations for subsonic wing-tail 
combinations . 


SYMBOLS 


A aspect ratio 

a velocity of sound in free stream . 

B = ^M 2 - 1 
b wing span 
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C I& = 


lift coefficient. 


lift- 

\ pV 2 Sw 


SCj 

3a 


3c 7 


CL ^ 0 


a ££ 

2V 

3Cr 


a — >0 


..-2 


a — >0 


horizontal-tail contribution to Cj- 


horizOntal-tail contribution to 






3c 7 
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2V 

3c T 


q — ^0 


a &£? 


q -*0 


C Iq = 


3c 7 


-3 

a S£_ 

8v3 


AC 




q .— >0 

horizontal-tail contribution to C 


Cm 


C n>a “ 


pitching-moment coefficient. 


Moment 

I pv 2 ^ 


^Cj^ 

3 a 


a — ^0 
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Sx " 
Cn£ = 

ACn^ 
C mq = 


^Cjh 

n, ac 
* 2V 

SCjh 

..-2 

■n ac 


la — >0 


a — >0 

horizontal-tail contribution to C^ 


horizontal-tail contribution to 


BCjh 


A qC 
° 2V 


q — >0 


^q " 


Sc 


-m 


% 

Cp 

c 

c 

d 

d o 

ds 


^Cm 


q — >0 


.—3 
B SSg. 
8v3 


q — >0 

horizontal-tail contribution to C^ 


Pressure 
\ pV 2 


pressure coefficient, 
root chord 

mean aerodynamic chord 

arbitrary distance along x-axis associated with a shift in 
center-of-gravity location (see table I) 

distance frcm. apex of a triangular wing to axis of pitch 
differential area 
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E 


E' 


complete elliptic integral of second kind with modulus k. 



complete elliptic integral of second kind with modulus k' 


F . force 


*0 


-A 


*2 


'0 


7 

G 

h 

h^_, ^2 


— times force associated with first term of equation (l) 
a 


— times force associated with second term of equation (l) 
ac 


4v 2 

..-2 

ac 


times force associated with third term of equation (l) 


— times force associated with first term of equation (2) 
qc 


4V 2 

. -2 
qc 

8y5 

—3 

qc^ 


times force associated with second term of equation (2) 


times force associated with third term of equation (2) 


arbitrary function associated with local angle-of -attack dis- 
tribution on a wing (see eq. (l4a)) 


finite part of an integral 


function associated with Mach surface formed by envelope of 
after Mach cones springing from trailing edge of wing 

arbitrary function associated with local angle of attack of 
an airfoil (see eq. (Dl)) 

limits of integration (see eq. (21)) 

function associated with equation of trailing edge of wing 
(see eq. (20)) 


J 
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K 


complete elliptic integral of first kind with modulus k, 
’' 1 dA 


X 


K q = -^ 


0 ]j( 1 - A 2 )(i - l^A 2 ) 

q 


rt 


L-_^3 2 - E* (Bn) + B2m - 2 . ■ K' (Bn) 

1 - B 2 * 2 1-bV 


K' complete elliptic integral of first kind with modulus k' 

k variable of integration 



*2 " 


x - x q 
Bnx q 


k 3 - 
kip = 


x - c 
Bnc 

Bmc 
x - c 



Z 

M 

Ms 

M 



distance from center of gravity to centroid of area of tail 

V 

free -stream Mach number, — 

a 

surface formed by envelope of after Mach cones springing from 
trailing edge of wing 

moment 

— times moment associated with first term of equation (l) 
a 

— times moment associated with second term of equation (l) 
ac 


D 


C 


-J 


(A 
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M2 


hi 2 - 

.. -2 

ac 


times moment associated with third term of equation (l) 



— times mcanent associated with first term of equation (2) 
qc 


— - times moment associated with second term of equation (2) 

• — d. 

qc 



m 


— - times moment associated with third term of equation (2) 

••-3 

qc 

tangent of semi vertex angle of triangular wing 
cotangent of trailing-edge sweep angle of wing 


P 

q 

q = 


• • 

q = 


Sq 

St 

Sfq 

St 2 


local static pressure minus free -stream static pressure 
rate of pitch 


R = fc 2 - B 2 (y 2 + z 2 ) 

S area of airfoil surface 

t time 

t Q arbitrary instant of time 




S9p 

discontinuity in — - which is induced within plan-form 

Sx 

boundaries by discontinuity in Qq across region of wake 


V free -stream velocity 

w velocity of vertical translation 
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X = x - xj 


Y = y - y-L 


x, y, z rectangular coordinates 

x l> Yl auxiliary rectangular coordinates 

x distance from leading edge of wing to center of gravity meas- 

ured in streamvd.se direction 

z distance between plane of wing and plane of tail 


a angle of attack, angle between body axes and free-stream 

direction at the center of gravity (see fig. 25(b)) 



the motion of center of gravity associated with this term corre- 
sponds to a constant vertical acceleration and is sometimes 
referred to as a plunging motion 


a = 


s? 


r spanwise circulation for an airfoil which has a spanwise 

circulation equal to ft(x,y) at wing trailing edge 

9 component of potential function resulting from constant 

vertical acceleration 


01, 0 2 


components of 0 


potential function 


potential function 
axis x = 


due 

due 


to constant rate of pitch 
to constant rate of pitch about the 


8 9 q 
8z qCy- 


average value of 


'N Q 

— at horizontal-tail location 

dz qc^. 


A variable of integration 

^ variable of integration 
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P 


a 



a l> a 2 


a 2 


0 

fo 


mass density of air 

local angle -of -attack distribution 


components of effective angle of attack on a tail which are 
associated with constant vertical acceleration 

average value of a 2 at horizontal-tail location 
potential function 

potential function associated with first term of equation (l) 



potential function associated with se.cond term of equation (l) 



potential function associated with third term of equation (l) 


4 

4 

4 




s 

^w 


potential function associated with first term of equation (2) 

potential function associated with second term of equation ( 2 ) 

potential function associated with third term of equation ( 2 ) 

components of potential function resulting from constant 
vertical acceleration (see eq. (9)) 

part of function ft associated with induced effect on wing 
part of function ft associated with wake 

potential function associated with constant angle of attack 


average value of — at horizontal-tail location 

dz aV 

Subscripts: 

LE leading edge of airfoil 


5 fyx, 

dz aV 
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TE trailing edge of airfoil 

w wing 

Expressions such as and 

derivatives of a tail surface considered as an isolated wing. Such 
expressions as u(o + ) mean that the value of SI is obtained by 
approaching zero from a positive direction. 



STATEMENT 'OF THE PROBLEM 


This paper contains an investigation of the contribution of the 
horizontal tail to the lift and pitching moment produced by a constant 
angle of attach, a constant rate of pitch, and a constant vertical 
acceleration. Unfortunately, the need for finding the force and moment 
associated with a constant vertical acceleration is not always recognized. 
Tn order to clarify the reasons for investigating this force and moment, 
an approach to the more general problem of calculating the effect of 
unsteady motions on the forces and moments on an aircraft undergoing a 
certain type of longitudinal motion is presented. 

Throughout the text the vertical velocity produced behind wings with 
a constant vertical acceleration at an angle of attack or in a pitching 
motion will be referred to as upwash. Negative upwash is, of course, 
commonly referred to as downwash. 

Consider the problem of finding the force and moment acting on the 
wing and horizontal tail of an aircraft undergoing longitudinal motions 
in which the forward velocity is constant. The force and moment are 
determined not only by the instantaneous velocities but also by their 
past values. The approach to the problem of calculating the effect of 
unsteady longitudinal motions presented in this paper takes into account 
the time history of the velocities as well as their instantaneous values. 

For many stability studies, the aircraft can be considered a rigid 
body and if the changes in the forward velocity are neglected, its longi- 
tudinal motion can be broken down into two components: (l) the vertical 

translation of the center of gravity and (2) the rotation of the center 
of gravity. 


The velocity of the vertical translation w and the velocity of 
rotation q of the center of gravity may be expressed as a power series 
in time about some arbitrary instant of time t Q . These two power series 
are 


| = a |t=t c + ( t 



(t - t Q ) 2 
21 


a 


t-t 0 




( 1 ) 
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and 


a(t) = q 


t=t r 


+ (t - t Q )q 


t=t r 


, (t - to) 2 •: 

i q 

21 


t=t r 


( 2 ) 


If it is assumed that the forces and moments can be calculated with 
sufficient accuracy by the linearized theory of supersonic flow, the 
following partial differential equation must be satisfied in the flow 
field: 





2V 

a^ 




tt 


= 0 


( 3 ) 


Since the preceding differential equation is linear, the flow associated 
with each term of equations (l) and (2) can be considered separately, 
and the total force and moment may then be found by summing the results 
from the individual terms. 


The boundary conditions associated with each of the first three 
terms of equations (l) and (2) are: Hie potential function associated 

with each term is zero upstream of the wing disturbance. On the wing and 
tail surfaces, the flow must be tangent to the surfaces; thus. 


dz 


-a,V 


Sz 


-qx 


c)z 


-dV(t - t Q ) 


¥1 

§ 7 " 


-qx(t - t Q ) 



dz 


aV 

21 


(t - t 0 ) 2 


2 _ _ qx / t 

Bz ” 21 


t o) 


2 
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In the region of the wake of the wing upstream of any disturbance from 
the tail surface, the pressure must he zero; thus, in the region of the 
wake 


< 1 _ n 
Sx v St 

300 + 1 30o 

Sx v St 

a £ + i a £_ 0 

301 + ! 30j 


Sx V St Sx V St 


+ 

Sx v St 


0 


Sx v. St 


0 


Solutions of equation ( 3 ) which satisfy the boundary conditions set 
forth in the preceding paragraph allow the calculation of the force and 
moment associated with each term of equations (l) and (2). The total 
force and moment can he written as 



/ - \ 2 


M 


= a lt=t 0 ^ + ° , |t=t 0 ^ *?(*) + ^|t=t 0 yy + 


+ 
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The F and M expressions in the preceding equations are functions 
of time and are not zero when t = t 0 . The force and moment at t = t 0 
can "be written in coefficient form as 


c X2 


C L = a |t=t 0 C Ia + a |t=t 0 2 V Cr h + a l t=t oUvj Cl u + * 


r\2 


fo\ 


q ^ § % + q C L^ + ^|t=t 0 ^ + * * * 


/ c \ 5 r 


00 


Cm = a |t=t 0 °nia + °'|t=t 0 °“a + a |t=t 0 ( 2 v) Si + * 


• • + 


q t=t n 2V 


S 3 . + q |-t^t 0 (|/) + q | t=t 0 ( 2 v} Cm q + 


(5) 


The coefficients C-^, C^, C^, and are the standard 

stability derivatives. The remaining quantities can he considered as 
stability derivatives which take into account the unsteady-flow effects 
other than the a effect. 

The series given by equations (4) and (5) are assumed to converge; 
however, in subsonic flow there is some evidence that the series may no 
converge for two-dimensional airfoils. In supersonic flow, the force 
and moment (calculated on the basis of linearized flow theory) depend 
only on the flow field between the leading edge of the wing and ^e 
trailing edge of the tail. Since only a small portion of the flow fiel 
is involved for supersonic flows, it is expected that these series will 

converge. 

In the present paper, it is assumed that the force and moment can 
be approximated sufficiently accurately by 



<xc 

2V 




( 6 ) 
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and 

^ = “S, + iSi + #S (7) 


The problem to he considered is the evaluation of the contribution of the 
horizontal tail to the derivatives contained in equations (6) and (7)» 


CALCULATION OF FLOW FIELDS OF AIRFOILS HAVING LINEAR 
ANGLE-OF-ATTACK VARIATIONS WITH TIME 
Basic Theory 


In the preceding section the aerodynamic lifting force and pitching 
moment were expressed as an infinite series of stability derivatives. 

The first stability derivatives of the series which are associated with 
unsteady flows were shown to be associated with the lift and moment 
resulting from a constant vertical acceleration. In this section methods 
are presented for the calculation of the flow fields produced by wings 
with a constant vertical acceleration and the potential function on the 
surface of tails which are located behind these wings. 

In order to calculate the flow field produced by a wing which has a 
constant vertical acceleration and the potential function on a horizontal 
tail behind the wing, the solution of two unsteady -flow problems will be 
utilized. The first problem is the calculation of the upwash induced 
by the wing at a constant vertical acceleration, and the second is the 
calculation of the pressure on an airfoil which has a local angle-of- 
attack distribution which varies linearly with time. These two problems 
axe special cases of the more general problem of calculating the flow 
induced by an airfoil which has an arbitrary angle -of -attack distribution 
which varies linearly with time. This general problem is formulated and 
the solutions to the other problems 'are taken from it. The effects of 
the rolling up of the wake and other distortions axe neglected and the 
wake is assumed to remain in the plane of the wing. 

The linearized partial differential equation (eq. ( 3 )) for unsteady 
supersonic flow is 

® 0xx " 0yy " 0zz + + 72 ^tt = ® 

S» EL 

where the x-axis lies along the free-stream direction. The potential for 
a zero -thickness lifting airfoil is antisymmetric with respect to the 



NACA TE 3072 


15 


plane of the wing (z = 0 plane). It is therefore necessary only to 
determine the potential above the plane for z = 0. This potential can 
be determined by specifying the boundary conditions in the plane of the 
wing and the potential upstream of the airfoil disturbance. The boundary 
and initial conditions for the potential for the region above an airfoil 
which has a local angle of attack varying linearly with time may be given 
as 


0=0 (8a) 

upstream of the wing disturbance, 

0 Z = -crVf (x,y) (t - t Q ) (8b) 

on the plan form, and 

P = -p(V0 x + 0t) = 0 (8c) 

for z = 0 (the plane of the wing) off the plan form. The condition 
given by equation (8c) is necessary to insure that there is no pres- 
sure discontinuity across the wake. 

From the preceding boundary condition, it follows that 
V0 x (x,y,O,t) + 0 t (x,y,O,t) = 0 

off the plan form. The general integral of this partial differential 
equation Is any arbitrary function of y and x - Vt. The potential 
is assumed to be continuous in the stream direction, and the potential 
has been taken to be zero upstream of the wing disturbance; thus, the 
function must be zero for points not on the plan form or in the region 
of the wake. The boundary conditions can now be written as 

0(x,y,z,t) = 0 

upstream of the wing disturbance, 

0(x,y,O,t) = 0 

for z = 0 for points not on the plan form nor in the region of the 
wake, 

0 z (x,y,O,t) = -crVf (x,y) (t - t D ) 
for z = 0 for points on the plan form, and 
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V0 X + ft = 0 


for z = 0 for points in the region of the wake. 

A solution of equation ( 3 ) can he found by assuming a function of 
the form 

0 = e(x,y,z) + xt(x,y,z) + ttt(x,y,z) ( 9 ) 

This assumption appears reasonable since the boundary conditions are 
linear functions of time and equation ( 3 ) has the same form if x 
and t are interchanged. 

Substituting equation ( 9 ) into equation ( 3 ) and equating powers 
of t yields 


B ^xx “ ^yy ~ ^zz ® 

B 0 xx + ^x + ® -^xx “ 0 yy “ '^yy - 9 zz “ x1 ^zz + p ^x “ ® 

If ijr is set equal to — 0, equation (10b) reduces to 

B 2 a 2 


(10a) 

(10b) 


® 20 xx “ 0 yy “ 0 zz — -z~ r(® 2fl xx - %y - n zz) - 0 


(H) 


Substituting equation (10a) into equation (ll) yields 


B 9 xx - 9 yy " e zz = 0 


( 12 ) 


Thus, under the assumption that 


\jr = — a 

B 2a 2 


equation ( 9 ) may be written as 

0 = 0 (x,y,z) + (t - ^ 0 n(x,y,z) 


( 13 ) 
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where 0 and Cl are solutions of the partial differential equation for 
steady supersonic flow. Note that Cl does not have the dimensions of 
a potential of steady flow. Equation (13) was found hy Clifford S. 
Gardner of New York University in an unpublished analysis using another 
approach and was used by Ribner in reference 16. 

The boundary conditions require that 

f z - e z + (^ - = “ t o) 

on the plan form or, if the powers of t are equated, 

Cl z = -oVf(x,y) (lka) 

and 

0 z " % = *Vt Q f (x,y) (14b) 

B^a 

Substituting equation (l4a) into equation (lkb) yields 

6 Z = (vt 0 - ^jrf(x,y) (Ike) 

The boundary conditions also require that 

0 = 0 (15a) 

Cl = 0 (15b) 

for all points either upstream of the wing disturbance or in the plane 
z = 0 which are not on the plan form or in the region of the wake and 
that 



Cl 


-■%" 0 


B 


ve x + 


vt 


(16) 
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for points in the region of the wake. When the powers of t are 
equated, equation (l6) yields 


= 0 


(17a) 


ve x - 



0 


(17b) 


Equations ( l4a) , (15b), and (17a) are the boundary conditions for 
a potential fl for an airfoil which has a local angle of attack equal 
to q f(x,y) in steady flow. Thus, since ft satisfies the linearized 
partial differential equation for steady flow, it can he determined hy 
the methods commonly used for obtaining steady-state solutions in super- 
sonic flow. 


Similarly, equations (l4c) and (15a) are the boundary conditions 
for the potential 0 for an airfoil which has a local angle of attack 


equal to 



xM^ cr f(x,y) 

B 2 J V * 


Also, because 0 satisfies the line- 


arized partial differential equation for steady flow, in the region 
unaffected by the wake it can be found by the methods used in steady- 
flow calculations. However, inasmuch as the boundary condition of equa- 
tion (17b) replaces the condition that in the region of the wake 0 X = 0, 
which is implied in these methods, another approach must be used to 
obtain solutions in the region affected by the wake. 


A convenient approach is to divide 0 into two parts: 0 = 0 B + Qq 
where 0p is the steady-flcw solution for the given airfoil in the 

ordinary sense, that is, subject to the condition that s- — =0 in the 

ox 

region of the wake. It follows that ©2 must satisfy the following 
boundary conditions: 

02( x ?y^ z ) = 0 

upstream of the wake influence, 

de 2 (x,y,z) 

— S = 0 

oz 


(18a) 


(18b) 


on the wing plan form, 
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0 2 (x,y,°) = 0 


(l8c) 


for points not on the wing plan form nor in the region of the wake, and 


de 2 (x,y,0) _ jL_ 

Sx B 2 V 

for points in the region of the wake. 


fl(x,y,o) 


(18a) 


Th e solution for 0 2 (x,y,z) can he obtained in several ways. Two 

8o 2 

approaches are presented here. Both methods require that ^ — he known 

ox 

- 

in the plane of the airfoil. The discontinuity in is zero every- 

where except in the region of the viscous wake, and except (for the case 
of wings with subsonic trailing edges) in the region of the wing affected 

8 0o 

by the wake, where a discontinuity in r-=- is induced by virtue of the 

ox 

conditions expressed in equations (l8b), (l8c), and (l8d) . The value of 


the discontinuity in 


80 r 

ST* 


for points in the region of the viscous wake 


can be determined from equation (l8d). In many cases, the discontinuity 
^0 

in induced on the plan form can be found by using the method devel- 

ox 

oped in reference (20) for the pressure cancellation for subsonic trailing 
edges. 

The first of the methods to be presented deals directly with the 


B0. 


in the z = 0 plane. If 


Z^0 2 


is assumed to be 


discontinuity in ^ ^ ^ — 

known in the plane of the wing, then 0 2 is given by (from eq. (6) , 
ref. 15 ) 


92( x >y> z ) = 


2*B 2 V 


K 


An(x 1 ,y 1 ,0)X 


Wake 


(y 2 + z 2 )r 


dXi dy-j^ + 


1 rr z An r X 

2rr JJ wing (y2 + z 2 )r 


dx 1 dy-L 


( 19 ) 
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Sgo 

where Au r is the discontinuity of - — across the z = 0 plane induced 

ox 

within the plan-form "boundaries by the discontinuity of 9p across the 
wake and where A ft is the known discontinuity in ft across the plane 
z = 0. Note that for wings with supersonic trailing edges Au r is zero. 

Since ft is the potential for an airfoil in steady flow, and since 
the potential discontinuity in the region of the wake for an airfoil in 
steady flow is constant in the x-direction, it follows that Aft is 
constant in the x-direction. In other words, Aft(x,y,0) is a function 
of y only and may be written as Aft(y). 

If the equation of the trail i ng edge is given by 

xq = J(yi) (20) 

equation (19) can be written as 


0 2 (x,y,z) = 


z r h 2 Aft(yq) nx-B^2+^2 

2kB 2 V J hq Y2 + z 2 dyi Jj( yi ) 



+ 


1 rr z Au r x 

s i/wing (y2 + z 2 ) e 4X1 dyi 


( 21 ) 


The limits on the outer integral of the first term of the right-hand 
side of equation (21) are taken across the span of the wing or to the 
limits of the forward Mach cone from the point (x,y,z). The first term 
on the right-hand side of equation (2l) can be integrated with respect 
to xq. The result of this operation is 


0 2^y> z ) = — ~2~ [ 2 ^ (yi) 2 i( x - J) 2 - B 2 (y 2 + z 2 )dy 1 
2*B 2 vJhq Y d + z 2 


q IT z AUj. X 
2* JJ Wing (y 2 + z 2 )r 1 yi 


( 22 ) 
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The second method of calculating 0 2 is similar to the approach 


used by Ribner in reference 1 6 . By this method an expression for 
is given and 0 2 is then found by integration. 


c^0 2 

dx 


The boundary conditions for 
trailing edges are 


d9g( x ,y,z) 

dx 


for -wings with supersonic 


de 2 (x,y_,z) 

dx 


upstream of the wake influence, 

d0 2 (x,y,O) 

Sx 


for points not in the region of the wake, and 


ae 2 (x,y,o) 

c)x 


- 5 - a( x ,y,o) 

B v 


for points in the region of the wake. These boundary conditions are 
the same as the boundary conditions for the potential of a lifting line 
located at the trailing edge of the airfoil. The strength of this 

lifting line is — r where P is the spanwise circulation for an 
B^V 

airfoil which has a spanwise circulation equal to the value of fl(x,y,z) 
at the wing trailing edge. Thus, 


d Q 2(x,y,z) 

dx 


1 

B 2 V 


(x,y,z) 


where ft w (x,y,z) is the part of ft(x,y, z) which results from the dis- 
continuity of ft(x,y,z) through the region of the wake. The function 
f2„.(x,y,z) is also the potential resulting from a lifting line which has 

a circulation of the amount P and is located at the trailing edge. It 
follows that 


1 r x 

e 2 (x,y,z) = — / dA 

B^V^Mg 


(23) 
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where the lower limit is taken from the trailing-edge Mach surface Ms. 
(The expression "trailing-edge Mach surface" refers to the surface formed 
hy the envelope of the after Mach cones springing from the trailing edge.) 

The solution for 0 2 (x,y,z) associated with wings which have sub- 
sonic trailing edges must include the effect of the potential induced on 
the plan-form area because of the condition that 


de 2 (x,y,o) 

dz 


on the wing plan form. Thus, if ftp denotes the effect of this induced 
potential for wings with subsonic trailing edges, 

0 2 (x,y,z) = -i- f fftp(A,y,z) + flw(A,y,z)1dA (24) 

B^V J Mg L - 1 

Two methods of finding the function 0 2 (x,y,z) have been presented, 
the expressions for 0 2 being given by equation (22) and equation (24). 
The function 0 2 (x,y,z) can also be obtained by other approaches, one 

of which will be utilized in an example which is presented at the end of 
this section. 

In appendix A, equation (22) is differentiated and the following 
equations are obtained: 


d9 2 (x,y,z) 

dy 


A^(yp) 


2jcB 2 V 


2(x - J) 2 


Y^ + z e 


+ B 


Jhp (y 2 + z 2 ) ]j(x - J) 2 - B 2 (y 2 + z 2 ) 


dyp 


z s rr x 

2rt g (y 2 + z 2 )r 


dxq dy x 


(25) 
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S9 2 (x,y,z) 

6z 


-vf 2 

2xB‘?V J h 1 


Afl(yi) 

(?* * zrf 




+ z 




z 2 (x - J) 2 

J) 2 - B 2 (y 2 + z 2 ) 


ayi + f - 1- 

2it oz 



X AUj. 
(y 2 + z 2 )r 


dxi <3yi 


(26) 


The differentiation of the double integrals of AUp was indicated only 


in the preceding equations , the reason being that the expressions 
resulting from the differentiation under the integral sign appear to be 
of little value in numerical calculations. For wings with supersonic 
trailing edges, AUj. is zero. Thus, the double integrals vanish. For 


most problems where the wing has subsonic trailing edges, the sidewash 
and upwash contributed by the double integrals in the usual tail loca- 
tion could be approximated by the sidewash and upwash from a lifting 
line in the vicinity of the trailing edge. This lifting line should 


have a strength equal to — i- (A0 2 ) ^ where (A0 2 ) ^ is the potential 


B 2 V 


'TE 


induced at the trailing edge by the discontinuity in 0 2 in the region 


of the wake. 


60p 

In appendix A the following expressions are developed for - — and 

oy 

— — for wings with supersonic trailing edges: 

6z 


i e g. (X -'Ai = 1 r « - 4- MMs) (27) 

6y B 2 V^G(y,z) dy B 2 V 6y 


a . 9g(X<y / Z) - 4- r ® - 4 - M*8> ( 28 ) 

6z B 2 V^G(y,z) ^z B 2 ? 6z 
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where and are the sidewash and upwash from a lifting line 

by bz 

with a strength equal to F and located at the trailing edge of the 
wing. For subsonic trailing edges. 


de 2 (x,y,z) 3_ P x dnp(A,y,z) 


dy 


~ f 

B 2 V 


dy 


dA + 


JLf* 

B 2 v 


x d%(>>>y> z ) 


dy 


dA (29) 


and 


c)z 


B0 2 (x,y,z) _ 1 T x Sflp(A,y,z) 

Jm_ Bz B 2 V^Mg dz 


B 2 V ^z 


d%(A,y,z) 


:/ 


dA (50) 


Bop bn t, 

where and — — are the sidewash and upwash resulting from the 

by bz 

potential induced on the plan form by the discontinuity of fl through 
the wake. 

The total potential (from eq. ( 13 ) ) in the region affected by the 
wake is, therefore, 

0(x,y,z,t) = 6 1 (x,y,z) + 0 2 (x,y,z) + cr^t - (31) 


The sidewash and the upwash are therefore given by 


y dy by \ B 2 a J by b 

and 


(32) 


B0^ B© 2 


K = t- + * - 


xM \ B ft 


z Bz Bz \ B 2 a ) bz a 


(33) 


The upwash on the plan form of a wing undergoing constant vertical 
acceleration is given by 
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0 Z = niV(t - t 0 ) 0 

for axes attached, to the wing. Therefore, for this motion, 

f(x,y) = 1 


a = a, 


The potential 0-^(x,y,z) is the steady -flow potential for an airfoil 
which has an angle -of -attack distribution equal to -^(vt Q - ^-jf(x,y) 

or ~lr^o “ and ^ or = 0 in the wake. Since a wing 

v \ B 2 / ° x 

having a constant rate of pitch about an axis located at x = x Q has an 

q(x - x n ) 

angle-of -attack distribution equal to — , 


9l(x,y,z) ^ 9q 


where 0^ is the steady-state potential due to a wing pitching about an 

B 2 Vt 0 

axis located at x = . 

M 2 

fl(x,y,z) 

The potential ; is the steady -flow potential resulting from 

a 

a wing having an angle of attack equal to f (x,y) or in this case unity. 

If is the potential due to a wing having a constant angle of attack a. 


«(x,y,z) = ^ % 

The function ©2 is given by equation (22) or equation (2k). 
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Therefore, the sidewash and upwash are given by 


and 


<h 


if d 9q 
B 2 by q_c 



to) 


M^x 6 2q | 1 692 
B 2 c by aV ac by 


(34) 


K 

ac 

where 0 q is 
Equation (35) 


if _S_ Sq 

B 2 6 z qc 




6 % ^ 1 d 0 2 

B 2 c bz aV ac 6 z 


(35) 


the potential for the wing pitching about the axis x = 0 . 
is used to calculate the upwash from a number of wings. 


Upwash at the Trailing Edge 

For wings with supersonic trailing edges, the upwash immediately 
behind the trailing edge is given by 


K 


if b_/0 + 

B 2 dzWw rpE; 


V 

C 


(t - t Q ) 


ifx 



I +i| 


| 

B 2 c 

dz' 

la vy 

CLC 

TE 

V'Jx J 

1 

TE 


(36) 


The first two terms of this expression may be evaluated by the method 
of calculating tbe steady-flow upwash immediately behind the trailing 

edge (refs. 10, 11, and 15 ) . The term ( ^-^) is (from appendix A) 
q \oz /q® 

the expression — — ^(Mg) — ^ Z - evaluated in the plane for z = 0 

B^ 6 z 

at the trailing edge of the wing. 


Approximate Values of 0y and 0 Z 

An approximation to the sidewash and upwash can be obtained by 
neglecting the terms which become constant for large values of x in 
equations (34) and (35)* This approximation leads to the following 
expressions: 
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h. « 



to) 


X % 

C dy aV 


(37a) 


& 

&c 


« 



to) 


x _S_ 

c dz aV 


(37b) 


Equations (37) niay also be obtained by assuming that the lift on the wing 
builds up instantaneously and that the upwash at a relatively large dis- 
tance behind the wing is determined by steady-flow values of the lift 
distribution on the wing. It is not assumed, however, that the upwash 
at the tail location at a given instant corresponds to the lift distri- 
bution on the wing at the same instant. The assumption is made that 
there is a time lag between the lift distribution on the wing and the 
upwash induced by this lift distribution at the tail location. This time 
lag is the time required for a point moving with the free -stream velocity 
to travel from the wing location to the tail location. In the past it 
has been found that, by using the stability derivative Cm^ calculated 

on the preceding assumptions, the airplane motions could generally be 
predicted with satisfactory accuracy (ref. 21). This accuracy seems to 
indicate that values of upwash from a wing having a constant vertical 
acceleration' in subsonic flow, calculated on the basis of the preceding 
assumptions, are a good approximation to the exact values. 


The question naturally arises as to whether the simple assumptions 
of . instantaneous lift build-up and time-lag effect give a good approxi- 
mation to the upwash from a wing which has a constant vertical accel- 
eration at supersonic speeds. This question is investigated for a number 
of cases. 


The values given by equations (37) are dependent upon the position 
of the origin of the x-axis. Since equations (37) represent a time-lag 
effect in the flow field, it seems logical that the origin of the x-axis 
should be located near the centroid' of the wing area for wings with small 
amounts of sweep. For highly swept wings it seems logical to allow the 
origin of the x-axis to vary along the span so that the distance x in 
equations (37) represents the distance from the point under consideration 
to the center of the section of the wing which lies directly upstream of 
the point under consideration. 

An approx im ate method of calculating the upwash behind wings with 
a constant acceleration which Is more accurate than equation (37b) is as 
follows: Calculate the angle-of -attach and pitching components approxi- 

mately by the vise of lifting lines. (See refs. 13 and 15.) The compo- 

S0p 

- — can be calculated exactly for wings with supersonic trailing 
oz 


nent 
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edges and appro xima tely for wings with subsonic trailing edges by neg- 
lecting the integral over the wing in equation (26). This method should 
yield results which are accurate to a fairly high degree except near the 
trailing edge or the Mach surface which springs from it since the main 
approximation is the calculation of the angle -of -attack and the pitching 
components by the use of lifting lines. 


The Two-Dimensional Airfoil 

One of the simplest time -dependent upwash problems is the two- 
dimensional airfoil which has a constant vertical acceleration. 

For the region which is not affected by the wake, 8g is zero, and 
by inspection 


80q . . 

- — = -q(x - Bz) 
oz 


The component 


8z 


is 


5% 

8z 


-aV 


From equation (35) , the upwash is given by 




(58) 


cS0q 8% 

From the region affected by the wake, — - and are zero and 

8z 8z 

equation (35) yields 


0 z _ JL_ 882 
dc dc 8z 


(59) 


From equation (28) , 


<V(Ms) L 


B^V 


8z 


8z 


(40) 
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When the y-axis is located along the leading edge of the wing (from 
eq. (A9) of appendix A) , 


G(y,z) = c + Bz 


m 


Since the potential function on the upper surface of a two-dimensional 

airfoil at a constant angle of attack is where x is measured 

B 

from the leading edge, ^(Mg) is given by 


<V(Ms) = 


q,Vc 

B 


( 42 ) 


Thus, the upwash in the region affected by the wake is found to be given 

■by 


0z = _l_ 
oc B^ 


m 


vz 

Values of — at t = t n are plotted for various values of z in 
a.c u 

figure 1. 

An examination of equation (37b) and figure 1 reveals that the 
approxi m ation given by equation (37b) breaks down for the upwash behind 
a two-dimensional airfoil. 


The Triangular Wing 

The upwash along the center line of the wake is determined for a 
triangular wing with subsonic leading edges. Equation (35) indicates 
that the upwash induced by a wing with a constant vertical acceleration 

is made up of three components. One of these three components, , 

oz_ aV 

is available in the published literature. (See refs. 11 and 12.) The 
two remaining components are derived in this paper. 

The upwash along the center line of the wake of a triangular wing 
with a constant angle of attack is given in references 11 and 12. From 
either of these references the upwash along the center line of the wake 
can be written as 
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JL 2 * 

Bz aV 


rtE'(Bm) 




E 


fc) 



+ 


r 1 K ( k i) - E fci) ^ + r 1 Kfe) - E(kg) 

J 0 k-L + Bm 1 ^^(l+k^) 



(Ha) 


for c ^ x ^ (Bm + l)c, and 


a Oa 2 

§z aV ” itE* (Bm) 


Ie K) 


^0 


K ( k i) - E ( k i) 


k]_ + Bm 


dk n 


(Hb) 


5 £Vr 

for (Bn + l)c ^ x. Figure 2 presents the variation of <r 2: along 

oz aV 

the wake center line of a triangular wing for various values of Bm. 
These values were taken from references 11 and 12. 

The upwash along the center line of the wake of a triangular wing 
with a constant rate of pitch is found in appendix B. Fran the results 
of this appendix the upwash along the wake center line can he expressed 
as follows: 


. _2K a c J 

dz 


~r 

c J 0 


ki[K(ki) - Efoj] 
(Bn + k 2 ) 2 


dkn 


: 


k 5 K(ia) - E(fc>) (l - k 5 2 )K(k 3 ) - E(kj) 

CjlKo - 

kg 2 ^Bmkg + l) 2 ^3 


(Ha) 
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for c < x < (Bm. + l)c, and 


de a (x,o,o) / x rH RbV + i)k x + E ( k ,) 


0 


(Bm + kq)‘ 


Bm + k]_ 


k'ikl + 


(2 - f)E0%) 


Em(l •- lq 4 . 2 )K(ki { .) 


k 4 


(45b) 


for (Bm + l)c ^ x. Figure 3 presents the variation of — ^ along 

qc dz 

the center line of the wake for various values of the parameter Bm. 

S8p 

The upwash component along the wake center line is deter - 

dz 

mined in appendix C. From the results of this appendix the upwash com- 
d9p 

ponent along the center line of the wake is 

dz 


dOg 

dz 


2%mc 


jtBE’(Bm) L 


2E(k 5 ) - (l - k 5 2 )K(k 5 ) 


(46a) 


for c S x =■ (ftn + l) c, and 


^02 _ 2a,(x - c) 

dz *B 2 E’ (Bm) 


2E(k 4 ) - (l - k4 2 )K(k4) 


(46b) 


for (Bm + l)c £ x. Figure 4 presents the variation of 
the wake center line for various values of the parameter 


b 2 d8g 

ac dz 
Bm. 


along 


The three components of the time -dependent upwash along the wake 
center line were added to yield the upwash at t = t Q . Figure 5 presents 

the results of this addition. The contribution of each component is 
shown in figure 6 for a Mach number of I .67 and an aspect ratio of 2. 
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Note that in this case no component is so small that it may he neglected. 
It can be seen from equation ( 35 ) that the upwash at any arbitrary time 
can be found from figures 2 and 5 * 

The results of exact and approximate calculations of the upwash 
behind triangular wings with a constant vertical acceleration are given 
in figures 7 and 8. Figure 7 presents values determined by the expres- 
sion (eq. (35)) for exact linearized upwash values and values determined 
by the approx imat e relation (eq. (37b)) for various aspect ratios for a 
Mach n umb er of 1.4l4. Figure 8 presents exact linearized values of the 
upwash and values determined by the approximate relation (eq. (37b) ) 
for various Mach numbers for an aspect ratio of 2. The origin of the 

2 c 

x-axis used in equation (3710 was located at the point. Figures 7 

and 8 indicate that equation ( 37 b) yields results which are good approxi- 
mations to the upwash along the center line of the wake of a delta wing 
with subsonic leading edges. 


The Rectangular Wing 


The upwash in the z = 0 p l ane behind a rectangular wing with a 
constant vertical acceleration is determined. From a knowledge of the 
upwash from the unswept wing of infinite span and the upwash from one 
tip of a semi -infinite rectangular wing, the upwash from a rectangular 
wing can be found as long as the Mach line from the leading edge of one 
tip section does not intersect the opposite tip. This was done for the 
wing at a constant angle of attack in reference 10 and for the wing 
pitching about its leading edge in reference 15 . 

The general upwash distribution behind a rectangular wing is not 
easily obtained in a convenient mathematical form. By use of equa- 
tion ( 36 ) the upwash close to the wing trailing edge can, however, be 
expressed in a concise mathematical form. The upwash close to the 
trailing edge of one tip of a rectangular wing is given by 


1 - 1 cos-l(aB + i)| 


•h 

dc 


itB £ 


cos- 1 /^ + 


sj.yfy *l)- 

rt c\c B / 


V(t - t 0 ) 


(Vf) 


where the coordinate axes Eire located at the leading edge of the tip 
section (see fig. 9)- The spanwise upwash distribution close behind the 
trailing edge of a semi -Infinite rectangular wing is obtained by putting 
negative values of y into equation (V7) . 

The upwash distribution for the rectangular wing is presented in 
the form of curves. Unfortunately, the results for the time -dependent 
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upwash from a wing with, a certain aspect ratio at a given Mach number 
can not be transformed to the results for the upwash from a wing with a 
different aspect ratio and at a different Mach number as can be done for 
wings in steady flows. The upwash for the motions considered here, how- 
ever, has been expressed in terms of steady -flow solutions. These steady- 
flow solutions can be transformed from one case to another. Each com- 
ponent of the upwash is given separately and then these components are 
combined for certain cases. 


The upwash from a rectangular wing at a constant angle of attack 
was given in reference 10. Figure 10 presents the upwash from one tip 
of a rectangular wing at a constant angle of attack. In figure 11 the 
upwash distributions from both tips are combined for various values of 
the aspect ratio Mach number parameter AB. 


The upwash from a rectangular pitching wing was presented in ref- 
erence 15. Figure 12 shows the upwash distribution from one tip of a 
rectangular wing pitching about its leading edge. In figure 13 the 
upwash distributions from both tips are combined for various values of 
aspect ratio Mach number parameter AB for the axis of pitch located 
at the wing midchord point. 


The upwash component 


302 

tr — was 
3z 


Figure 14 shows the distribution of 
one tip of a rectangular wing. 


evaluated by using equation (28) . 

•n2 30p 

the upwash component — - behind 

oc oz 


The three components of the time -dependent upwash for one tip of a 
rectangular wing were added to yield the time -dependent upwash behind 
one tip of a rectangular wing at t = t Q . Figure 15 presents the results 
of this addition for various Mach numbers. The contribution of each 
upwash component is shown in figure 16. This figure indicates that no 
one of the three components which make up the upwash at t = t D is so 

small that it may be neglected in upwash calculations. In figures 17 
to 22 the upwash distributions from both tips are combined for various 
Mach numbers and aspect ratios. From equation (35) it can be easily 
seen that the upwash at any time can be determined from figures 11 and 17 
to 22. 


Figure 23 presents exact (linearized) and approximate values of 
upwash at t = t 0 for a wing with an aspect ratio of 1 for Mach numbers 
of 1.4l4 and 2.24 for x/c of 3 and 6. Figure 24 presents exact (line- 
arized) and approximate values of upwash at t = t Q for a wing with an 
aspect ratio of 3 tor Mach numbers of 1.4l4 and 2.24 for x/c of 3 and 6. 
These two figures indicate that the approximate values are in fair agree- 
ment with the exact values for x/c of 6. As pointed out previously, the 
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flow behind a rectangular airfoil may be considered as being made up of 
the effect of the tips and the effect of a two- dime nsional airfoil. 

Since the approximation breaks down for the two-dimensional airfoil, 
the effect of the tips must be considerably larger than the effect of 
the two-dimensional airfoil in order that the approximation be valid. 

This explains the large discrepancy between the exact and approximate 
values shown in figure 24 for the region near the center of the wake. 

It should be noted that the upwash distributions shown in fig- 
tires 11, 13, and 17 to 2k have been plotted from the wing axis of symmetry 
out to the wing tip. In figures 10, 12, 14, 15, and 16, showing upwash 
distributions from one tip of a rectangular wing, the origin is located 
at the wing tip. 


AERODYNAMIC COEFFICIENTS 


In order to calculate the contribution of the horizontal tail to 
the aerodynamic coefficients a knowledge of the upwash produced by the 
wing is required. Exact (linearized) and approximate methods for the 
calculation of the upwash induced by a wing which has a constant vertical 
acceleration were presented In the preceding sections. It was found that 
the exact (linearized) upwash induced behind wings with a constant accel- 
eration contains the upwash induced behind the same wing at a constant 
angle of attack and with a constant rate of pitch. In this section exact 
and approximate methods for the calculation of the aerodynamic coeffi- 
cients on the horizontal tail are presented. 

The calculation of the tail contribution to the derivatives con- 
sidered in this paper can be accomplished by utilizing the concept of 
the "effective angle of attack. " The effective angle of attack is 
defined in such a way that the local single -of -attack distribution on the 
tail surface is changed in order to take into account the upwash induced 
by the wing. The tail surface is then considered as an isolated lifting 
surface on which the local angle -of -attack distribution is given by the 
effective angle of attack. Except for the two-dimensional wing-tail 
combination treated, the exact (linearized) values of the aerodynamic 
coefficients to be presented were obtained by numerical integrations. 

The coordinate axes used in determining the aerodynamic coefficients 
are illustrated in figure 25(a). The origin is located at the center of 
gravity of the aircraft. The effect of changes in the center -of -gravity 
locations can be taken into account by the transformations presented in 
table I. The stability axes are illustrated in figure 25(b), and the 
transformations from body axes to stability axes are given in table I. 
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The Lift and Moment Due to Angle of Attack for 
Wing-Tail Combinations 


Equations for tail contribution . - The effective angle of attack used 
in calculating the tail contribution to Ci^ and Cn^ is 


a 


a^l + 


AM 

8z aV / 


m 


The lifting pressure associated with equation (48) can be found by 
steady-flow methods. 

The normal-force and moment coefficients resulting frcm the pres- 
sure difference through the tail surface can be expressed as 


ACr = — fT AP ds (49) 

apV 2 S t Sw JJrpaii 


AC™ ff x AP ds 

apV^Stet ^Taii 


where x is measured from the axis about which the moment is taken. 

An approximation of the contribution of the tail to the stability 
derivatives and can be found by making use of an arithmetic 

average of the upwash induced by the wing at the tail location. When 
this is done, AC^ AC^ W be expressed as 


^Ia 




(51) 



!_ aY_l ct 1 

8z aV y\c w J ixlt 



( 52 ) 
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7 . 

For most cases -±- Cj^J^ is much greater than Cm^ | and there- 

c w °w 

fore equation (52) can he approximated by 





+ 


-6_ c T I 
6z aV/5^ ^It 


(53) 


Two-dimensional wing -tail combinations .- The and the C m ^ 

for the two- dim ensional wing-tail combination can be expressed in closed 
form. The sum of the wing and the tail contributions to these deriva- 
tives is as follows: For the case where the tail lies downstream of 
the Mach sheet from the trailing edge of the wing (see fig. 26(a)), 


and 





( 5*0 


(55) 


For the case where the Mach sheet from the trailing edge of the wing 
intersects the tail surface (see fig. 26(b)), 


and. 



(56) 








2xBz 



(57) 


For the case where the tail surface lies between the Mach sheet from the 
leading edge of the wing and the Mach sheet from the trailing edge of 
the wing (see fig. 26 (c)), 
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and 




(58) 


(59) 


Figure 27 presents a comparison of values of the calculated 

by the preceding exact expressions with the values of the Cm^ where 

the tail contribution was calculated by the approximate relation given 
by equation (53) • Equation (53) is a good approximation to the tail 
contribution to the for the case considered in figure 27. The 

discontinuities in the curves in figure 27 correspond to points where 
the type of flow over the tail surface changes. (See fig. 26.) 

A triangular wing-tail combination .- The and of a tri- 

angular wing-tail combination are considered. The wing has an aspect 
ratio of 2.31 and the aspect ratio of the tail is twice that of the 
wing. The wing and the tail are located in the same plane, and the root 
chord of the wing is four times the root chord of the tail. The Mach 
number is restricted to the range' where the leading edge of the wing 
is subsonic and the leading edge of the tail is supersonic. 

In reference 11 calculated values of upwash behind triangular wings 
with subsonic leading edges are presented. These values could be used 
to calculate the tail contribution to the C]^ and of the wing- 

tail combination which is being considered here. In order to be con- 
sistent with what follows, however, it is assumed that the spanwise 
variation of the upwash can be neglected and the upwash is determined 
by the values of the upwash along the center line of the wake. (An 
examination of the upwash values presented in reference 11 indicates 
that spanwise variation of the upwash is small in the region where the 
tail is located. ) 

Figure 28 presents the variation of with 

Mach numbers for the triangular wing-tail combination under consider- 
ation. Figure 29 presents the variation of with l/c^ for two 

of the Mach numbers considered in figure 28. Figure 30 presents a com- 
parison of the exact (linearized) values of the with the values 

of On ^ where the tail contribution was calculated by the approximate 
relation given by equation (53). Figure 30 indicates that equation (53) 


l/Cw 1 1 ° r three 
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is a good approximation to the tail contribution to for the cases 

being considered in this section. 

Rectangular -wing — triangular -tail combinations . - The and 

of a number of rectang ula r -wing — triangular-tail combinations are con- 
sidered. These wing-tail combinations are illustrated in figure J>1 where 
the defining parameters are also presented. The center of gravity is 
located at the midchord point of the wing. 

The variation of Cp^ with 2/% of rectangular -wing — 

triangular -tail combination for several Mach numbers is presented in 
figure 32 whereas the variation of Cm is given in figure 33* The 

decrease in C-^ with 

the upwash with increasing l/c ■ w . In certain cases, the effect of the 
upwash from, the wing on the tail is large enough to cause the total 
Cm-, to become destabilizing (fig. 33)* 

Figure 3^ presents a comparison of the exact (linearized) values 
of Cjj^ with the values of where the tail contribution was cal- 

culated by the approximate relation given by equation (53)- Figure Jh 
indicates that equation (53) is a good approximation to the tail con- 
tribution to for the cases considered. 


The Lift and Moment Due to Steady Pitching 
for Wing-Tail Combinations 

General expressions for the aerodynamic coefficients .- The effec- 
tive angle of attack of a tail surface in pitching motion is 


increasing 1/Cy is a result of the increase in 


^w/x- + _3_ 9 q. \ 

V 8z qcj 


(60) 


where the origin of the x-axis is located at the axis of pitch. The 
lifting pressure associated with a can be found by steady -flow methods. 

The normal-force and moment coefficients resulting from the lifting 
pressure on the tail surface can be expressed as 
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& = i ft rr £S? ds 

^ Sy JJ sur face 

2V 


(61) 


1 stct np 


Tail surface 


x ACp 

qpy 

2V 


ds 


(62) 


where the moment axis is taken about the axis of pitch. 


Approximations of the contribution of the tail to the stability 

C can be found by making use of an arithmetic 

mq 


derivatives and 


average of the upwash induced by the wing at the tail location. When 

— n 1 1 f J J 1 A /vrnU ■? -r -1 o4"^ r\Vl TT1 Q Tf h A 

this is done 
expressed as 


this is done, the and of the wing-tail combination may be 


C Lq = C Lq|wing + S 2 ^ ^Tail + 2 L qSy + c* %kil 


( 63 ) 


St 


% %|wing + 


_2 (t) ^i^ 11 ' h C ^k 11 


_L iL 9 ( 3 - 

Cy 3z qCy 


C *tjTail + 


2 


ft JL fl_ n I + 
Cy- dz qiCy ^“a-ITail 



Tail 


+ 


Cy S^Tail 


( 64 ) 


where in this case Z is 
centroid of the tail area, 
the stability derivatives 
of the tail area.) 


the distance from the axis of pitch to the 
(The axis of pitch and the moment axis for 
of the isolated tail are located at the centroid 
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Further approximations to Cj^ and Cj^ can he obtained by 

retaining only the highest powers of Z/Cy in the preceding expressions. 
The following approximate relations are thereby obtained: 


Cr « 2 — J- Or I 

^ SyCy ^1 


Tail 


(65) 


s 


« -2 


SwV 


Sf* I 


Tail 


(66) 


Two-d~i Tnpnsional wing -tail combinations . - The equations for C-^ 

nnri for the two-dimensional wing-tail combinations can be expressed 

in closed form. The sum of the wing and the tail contributions to these 
derivatives are as follows: For the case where the tail lies downstream 

of the Mach sheet from the trailing edge of the wing (see fig. 26 (a)), 


C I= |(1 - 2 — + 

H B\ Cy 


lc t \ 


-v 


(67) 


and 


% B 


°w 


r L) * i/lr + ip-V 

\ c w/ 1 ^V c w/ Cwvwv 


( 68 ) 


For the case where the Mach sheet from the trailing edge of the wing 
intersects the tall surface (see fig. 26(b)), 


Ct = — 
H B 


and 


5BZc t _ 2BZZ + -l _ 2x + |ft . + f mf _ (js^\ 2 + /j_\ 2 + / c t 

Cy 2 Cy 2 °w Cy 2 \ c w J \ c w j \ c wy ' 2C 1 


-W; 


(69) 


% " " B 


c t 


2 + x 2 


c w\3 


+ JL _ 1 + I 1 / 


5lc t . B 2 ^ 2 


-v 


Cv 


^W 2 5 c w 2 c w 2 


-i + a- 

°w 


^°w °w 


X 2 \ _Z_/k lf_ + + Ct2 ' 

2] + Cy\3 - 2 + . 2 + _ 2, 


ct- 


°w 


Cw 2c w 


12V 


( 70 ) 
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For the case where the tail surface lies between the Mach sheet from 
the leading edge of the wing and the Mach sheet from the trailing edge 
of the wing (see fig. 26(c)), 


and 




2x + 2Bzc t \ 
cv 2 / 



JL + /JL \ 2 + BzZc t 

°w VV/ 


(71) 


(72) 


Figure 35 presents a comparison between exact (linearized) values 
of the Cjq calculated from equations (67), (69), and (71) and approxi- 
mate values calculated by the use of equation (65) . Figure 36 presents 
a comparison between exact (linearized) values of the calculated 

from equations (68), (70) , and (72) and approximate values calculated by 
the use of equation (66). Figures 35 and 36 indicate that the approxi- 
mate relations given by equations (65) and (66) yield results which are 
in good agreement with the exact (linearized) values for if c^. greater 

than 2 for the cases considered. The poor agreement between exact and 
approximate values of Cj^ for the values of lj c^ less than I.75 

(fig. 35) are a result of neglecting the effect of the upwash induced 
by the wing on the tail surface . 

A triangular wing-tail combination .- The and of a 

triangular wing- tail combination are considered for the same configu- 
rations and Mach number range which were considered previously in the 
section on the lift and pitching moment due to angle of attack. 

It is assumed that the spanwise Variation of the upwash due to 
pitching can be neglected and that the upwash at the tail locations can 
be determined from the values of upwash along the center line of the wake. 
Under these assumptions the upwash at the tail location due to the wing 
pitching about its apex can be determined from figure 3. The upwash at 
any point on the center line of the wake due to the wing pitching about 
an axis at an arbitrary location can be found by the use of figures 2 
and 3 and the relation 


_ _1_ d 0( l 

Pitching about < l c 
arbitrary axis 


Aq 6 flq 
c ’bz aV 


JL 

qc 8z 


Pitching 
about apex 
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where do is the distance from the apex of the wing to the location of 
the arbitrary axis of pitch. 

The variation of and with for three Mach numbers 

is presented in figures 37 and respectively. 

A comparison between exact (linearized) values of C-^ and the 

approximate values of Cj^ calculated from equation ( 65 ) is given in 

figure 39 - The large discrepancy between the exact and approximate 
values for the lower Mach number. is a result of neglecting the upwash 
from the wing due to pitching and of the wing. 

Figure 40 presents a comparison between exact (linearized) values 
of the and approximate values of the calculated from equa- 

tion ( 66 ). Thi s figure indicates that, for values of l/c w greater 
than 2 , equation ( 66 ) yields a good approximation to for the two 

cases considered. 


Plots of Ciq and for three different center -of -gravity loca- 

tions at Mach number 1.44 are presented in figures 4l and 42, respec- 
tively; these plots indicate that the center-of -gravity location has a 
large effect on Cj^ and for the cases considered. 

Figure 43 presents a comparison between exact (linearized) values 
of and values calculated from equation ( 65 ) for two different 

center -of -gravity locations. The poor agreement between exact and 
approximate values shown in figure 43 indicates that equation ( 65 ) should 
not be used when the center of gravity is not located near the centroid 
of the wing. 

For two different center -of -gravity locations, a comparison between 
exact (linearized) values of and values calculated from equation (66) 

is made in figure 44. Since the agreement for the lower values of Z/c^ 

is poor, this comparison indicates that care should be exercised when 
equation (66) is used for cases where the center of gravity is not located 
near the centroid of the wing. 


Rectangular -wing — triangular -tail combinations . - The 


% ^ % 


of a number of rectangular -wing — triangular -tail combinations are con- 
sidered. These wing-tail combinations are illustrated in figure 31 where 
the defining parameters are also presented. These are the same configu- 
rations that were considered previously in the sections on the lift and 
pitching moment due to angle of attack. 
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The variation of Cp^ with if of the wing-tail combination for 
several Mach numbers is presented in figure 45 whereas the variation of 
Cmg_ Z/Otf 


is presented in figure 46. 


A comparison of the exact (linearized) values of 



and approxi- 


mate values calculated from equation ( 65 ) for two wing-tail configura- 
tions and two Mach numbers is presented in figure 47 . There is good 
agreement between the exact and the approximate values for the cases 
considered. 


Good agreem e nt between the exact (linearized) values of C^ and 

the approximate values calculated from equation ( 66 ) for the cases con- 
sidered is indicated in figure 48 which presents these values for two 
wing-tail configurations and two Mach numbers. 


Figures 49 and 50 present Cp^ and C^, respectively, for three 

different center -of -gravity locations at Mach number 1.4l4 for two dif- 
ferent wing-tail combinations. These figures indicate that the center - 


of -gravity location has a large effect on 
considered. 




and 


Cmq 


for the cases 


A comparison between exact (linearized) values of Cp^ and values 

calculated from equation ( 65 ) for two different center -of -gravity loca- 
tions and two different configurations is presented in figure 51» The 
poor agreement between the exact and approximate values shown in fig- 
ure 51 indicates that equation ( 65 ) should not be used when the center 
of gravity is not located near the centroid of the wing. 

Figure 52 presents a comparison between exact (linearized) values 
of the and values calculated from equation ( 66 ) for two different 

center-of -gravity locations and two different configurations. The poor 
agreement for the lower values of l fc^ shown in figure 52 indicates 

that results calculated from the approximate equation ( 66 ) are unreliable 
for center-of -gravity locations that are not near the centroid of the 
wing area. 


The Lift and Moment Resulting From a Constant Vertical 

Acceleration for Wing-Tail Combinations 

General expressions for aerodynamic coefficients .- The upwash on the 
horizontal-tail surface resulting from the constant vertical acceleration 
is given by 


0 Z = _aV(t - t Q ) 
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The upwash induced, "by the wing at the tail location is (from eq. (35)) 


a . M 2 d 8q 

= ac — ^ + ac 

z g2 az qc 


2 (t - t o) - 

c B 2 c 


d 2a ^ e 2 
dz aV §z~ 


The effective angle of attack of the tail is given "by 


. , + \(i ^ Pa\ ac M 2 d 9q ac ^ x S % 1 ^ 6 2 

° a 0 \ + dz aV/ + V p2 Bz qc V g2 c Sz aV V Bz 


The preceding equation indicates that the effective angle of attack of 
the tail is made up of two components, one which varies linearly with 
time and one which is of a steady-state nature. It follows that, in 
order to determine the potential induced by the tail on its surface, 
it is necessary to solve two problems, one which has an angle -of -attack 
variation given by 



(73) 


and one which has an angle -of -attack variation given by 

_ ac M 2 _d_ £q ac x d ^a + 1 ^02 
02 V p2 Bz qc V j2 c 3z aV V dz 


(74) 


The potential corresponding to cr^ can be found by the methods given in 

the theoretical development concerning upwash and the potential corre- 
sponding to 02 can 4> e found by steady -flow methods. 

After the potential is determined, the pressure can be found by 
making use of the relation 

AP = p(V Apx + Apt) 


where Ap represents the potential difference induced by the tail through 
its surface. 
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The normal force and moment resulting from the pressure difference 
through the tail surface can he expressed by 




AP ds 


Tail surface 


and 


-a 

M l = " 


x AP ds 


Tail surface 


The two quantities , and can be expressed in coefficient form by 




£ 


£ V^Sy pV 2 ^ ^Tail surface 


AP ds 


and 


-a 

M 1 


£ V^SyC^ p^SyCy JJla.il surface 




x AP ds 


The tail contribution to the stability derivatives 
be expressed as 


Ct, and Cm. 
"u, a 


can 


^14 - 


St^t rr 

apVS^c-t SyCy JJ Qjail surface 


AP ds 


( 75 ) 


and 




k s t^t 2 
apVStC t 2 Sy5y' 


M 


x AP ds 


Tail surface 


( 76 ) 


The difficulties involved in calculating the exact values of the 


tali' contribution to the stability derivatives C-^ and 




are in 
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most cases too great for practical computations. The lift on tails with 
supersonic leading edges and trailing edges which are perpendicular to 
the free -stream direction can, however, he treated hy the method pre- 
sented in reference 22. Expressions for the moment on the tail are 
developed in appendix D for the case of an angle -of -attack distribution 
that increases linearly with tine. 


An approximation to the contribution of the tail to the stability 
derivatives and Cm^ can be found (as was done in ref. l 6 ) by 

mnlcl ng use of an arithmetic average value of the upwash induced by the 
wing at the tail location. 

expressed by 


When this is done, AC-^ anl 


AT - StCt 
^ Sw c w 


( x + Sz hv) Cl d 


Tail 


,2V - - r I 
+ 5^ °2 C IalTail 


(77) 


and 





I + £1V_ — n I . 

C t \ Bz a V/ 1 Tail ac^ ^ At, I Tail 


1 + 




Tail 


2 Va p 

+ 


ac 4 


Cm 


a I Tail 


(78) 


where l is the distance from the center of gravity to the reference 
axis about which and 8X6 taken. 


Further approximations to 


C I& and 


can be obtained by 


retaining only the terms which are of the highest order in l in equa- 
tions ( 77 ) and ( 78 ). When these terms are retained, the following 
approximate relations are obtained: 



^t l_ 5 2q 
Sy- c^ 57 dV 


Cl tJTail 


(79) 


and 

^ “ 2 I I 


(80) 
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Equations ( 79 ) and ( 80) can also "be obtained by assuming a time-lag 
effect in the downwash from the wing. Approximate expressions for wing- 
tail combinations are 


and 


“ c ia.l 


Wing 


^t l_ £a £ j 

Sy.Cy.SzaV ^tolTail 


(81) 


Sx~ 2 


s t /] \ 2 S fig 
^w\ 9 w/ 8z aV 


\\ 


Tail 


(82) 


The of the wing has been neglected in equation (82) because it is 

usually small. 


Two-dimensional wing -tail combinations . - The equations for 



and for two-dimensional wing-tail combinations can be expressed 

in closed form. The sum of the wing and the tail contributions are as 
follows (these expressions were obtained by integrations of the pressure 
distribution resulting from a constant vertical acceleration over the 
wing and tail surfaces): For the case where the tail lies downstream 

of the Mach sheet from the trailing edge of the wing (see fig. 26(a)), 



For the case where the Mach sheet from the trailing edge of the wing 
intersects the tail surface (see fig. 26(b)), 


0 - _ JL 

b 3 


1 + 


, x_ _ -,\ (l_ _ Bz 

2 lCy Cy 2Cy + Cy ^ J (Cy Cy 


+ + Jbs. + J 2 ^ SBa t 

2Cy Cy Cy / n_2 
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(85) 
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For "the case where the tail surface lies between the Mach sheet from the 
leading edge of the wing and the Mach sheet from the trailing edge of 
the wing (see fig. 2.6(c)), 



(87) 





( 88 ) 


A comparison between exact (linearized) values of the and 

values calculated from equation (8l) is presented in figure 53 for two 
Mach numbers. The poor agreement between exact and approximate values 
shown in figure 53 is to be expected because the upwash behind the 
two-dimensional airfoil can not be approximated by a simple time-lag 
effect. 

Figure 54 presents a comparison between exact (linearized) values 
of the and values calculated frcm equation (82) for two Mach num- 

bers. This figure indicates that equation (82) can not be used for two- 
dimensional wing-tail combinations. The poor agreement between exact 
and approximate values shown in figure 54 is to be expected because the 
upwash behind the wing does not exhibit a time-lag effect. 
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A triangular -wing-tail combination .- The 


and 



of a 


triangular wing-tail combination are considered for the same configura- 
tion with the same Mach number range which was considered previously in 
the sections on the lift and pitching moment due to angle of attack and 
in the sections on the lift and pitching moment due to pitching. 


It is assumed that the spanwise variation of the upwash due to a 
constant vertical acceleration can be neglected and that the upwash at 
the tail location can be determined from the values of the upwash along 
the center line of the wake. Under these assumptions the upwash at the 
tail location due to a constant vertical acceleration at t = t Q can be 
found from figure 5- 


The variation of C]^ and with l /c^ for three Mach numbers 

is presented in figures 55 and 56, respectively. 

A comparison between exact (neglecting spanwise variations in 
upwash) values of and approximate values calculated from equa- 

tion (8l) (fig. 57) shows good agreement between the exact and approxi- 
mate values. 


A similar comparison between exact (neglecting spanwise variations 
in upwash) values of and approximate values calculated from equa- 
tion (82) (fig. 58) also shows good agreement between the exact and the. 
approximate values. 


Figure 59 presents the variation of for three different center - 

of -gravity locations for Mach number 1.4-4- . This figure indicates that 
the center-of -gravity location has a large effect on C^. 


Exact (neglecting spanwise variations in upwash) values of 


Si 


and values calculated from equation (82) for two different center-of- 
gravity locations at Mach number 1.44- are compared in figure 60. The 
agreement between exact and approximate values is much poorer in these 
cases than for the cases shown in figure 5^. This result indicates that 
equation (82) should not be used for cases where the center of gravity 
is not located near the centroid of the wing area. 


Rectangular -wing — triangular -tail combinations . - The and 

of a number of rectangular -wing — triangular-tail combinations are con- 
sidered. These wing-tail configurations are illustrated in figure 51 
where the defining parameters are also presented. These are the same 
configurations that were considered previously in the sections on the 
lift and moment due to angle of attack and in the sections on the lift 
and moment due to pitching. 
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The variation of Cj^ and with l j c^. for several Mach num- 

bers is presented in figures 6l and 62, respectively. 

Figure 63 presents a comparison "between exact values of and 

values calculated from equation (8l) for two different wing -tail combi- 
nations. Figure 63(a) shows good agreement between exact and approxi- 
mate values for a wing with an aspect ratio of 1 whereas figure 63(b) 
shows poor agreement between the exact and approximate values for a wing 
with an aspect ratio of 4 . This poor agreement is to be expected because 
the upwash behind the wing with an aspect ratio of 4 (fig. 63(b)) exhibits 
some of the two-dimensional characteristics. 


A comparison between exact values of C^ and the values calculated 

from equation (82) for two different wing-tail combinations is presented 
in figure 64 -. Figure 64 (a) (aspect-ratio -1 wing) shows fairly good 
agreement between exact and approximate values, and figure 64 (b) (aspect- 
ratio -4 wing) shows poor agreement which is caused by the two-dimensional 
character of the upwash behind the wing. 


The variation of C ^. with 2 /Cy for three center -of -gravity loca- 


tions for two configurations is given in figure 65. This figure indicates 

that the center -of -gravity location has a strong effect on Cm, . 

a 


Figure 66 presents a comparison between exact values of Cm. and 

approximate values calculated from equation (82). Figure 66(a) indi- 
cates that equation (82) should not be used for cases where the center 
of gravity is not located near the centroid of the wing area for a wing 
of low aspect ratio whereas figure 66(b) indicates that equation (82) 
should not be used at ail for wings of high aspect ratio. 


The Cm^ + Cm^ for a Number of Wing-Tail Combinations 

The expression C^ + C^ partly determines the damping of longi- 
tudinal oscillations of aircraft. For this reason, C^ + C^ is given 
separate consideration in this section. 

Figure 67 presents C^ + C^ for a series of two-dimensional wing- 

tail combinations in which the center of gravity is located at the mid- 
chord point of the wing, lj c^ is equal to 2.25, the ratio is 

i, and the height of the tail surface above the wing has various 
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values. The discontinuities in slope of the curves in figure 67 corre- 
spond to points where the type of flow over the tail surface changes 
(see fig. 26 ). Figure 67 indicates that, when the tail lies between 
the Mach sheets from the leading and trailing edges of the wing (see 
fig. 26(c)), + Cm^ is decreased considerably from the values for 

the same combination where the wing and the tail lie in the same plane. 


The variation of C^ + C^ with Z/% is presented in figure 68 

for two Mach numbers for a two-dimensional wing-tail combination in 
which the center of gravity is located at the midchord point of the wing, 
the ratio c-fc/cy- is i, and z/c w is K This figure indicates that 
I s decreased by increasing Z/c w . 


Figure 69 presents the variation of 


Cm^ + with Mach number 


for a two-dimensional wing-tail combination in which the center of 
gravity is located at the midchord point of the wing, the ratio l jc^ 

has various values, the ratio c^/e^. is i, and the tail lies in the plane 
of the wing. 


The relation between Mach number and center -of -gravity location 
which causes C^ + to be zero for five two-dimensional wing-tail 

combinations is shown in figure 70. 

Figure 71 presents the variation of C^ + C^ with Z/c w for 

three Mach numbers for the same triangular wing-tail combination as was 
considered previously. Figure 72 presents the variation of 

with Z / Cy for various Mach numbers for the same rectangular -wing- 

triangular -tail combinations as were considered previously. These fig- 
ures indicate that Ch + C^. increases rapidly with Z/c^. 


CONCLUDING REMARKS 


The force and moment coefficients of an aircraft undergoing unsteady 
longitudinal motions at supersonic speeds can be expressed in terms of an 
infinite series of stability derivatives of successively higher orders. 
This representation of the aerodynamic forces and moments is felt to be 
useful in accounting for the unsteady influences in stability studies. 

In this paper attention is primarily devoted to the establishment of the 
more common stability derivatives. The stability derivatives resulting 
from unsteady motions that appear to be' the most important are those 
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associated with constant vertical acceleration. The calculation of the 
upwash by the linearized theory behind wings with constant vertical accel- 
erations can be reduced to solving a number of steady -flow problems. 

One of these steady-flow problems is the determination of the upwash 
behind a wing at a constant angle of attack and another of these steady- 
flow problems is the determination of the flow behind a wing indu ced by 
a constant rate of pitch. The effects of the rolling up of the wake and 
other distortions are neglected and it is assumed that the wake remains 
in the plane of the wing. 

For some plan forms the upwash at points which are not located near 
the Mach sheet from the trailing edge of the wing resulting from a con- 
stant vertical acceleration can be approximated very well by using the 
upwash due to a constant angle of attack and a time-lag effect. This 
approximation breaks down for the two-dimensional airfoil and yields 
poor results for high-aspect-ratio rectangular wings. The few calcu- 
lated examples seem to indicate that in general a simple time-lag effect 
yields good approximations to the upwash due to a constant vertical 
acceleration behind unswept wings of low aspect ratio and yields a poor 
approximation to the upwash behind unswept wings of high aspect ratio. 

The effects of sweep were not investigated in any of the examples which 
were calculated. 

For the wing-tail combinations investigated, the results indicate 
that the moment coefficient resulting from a steady pitching can be 

approximated to a fairly high degree of accuracy by a simple expression. 
This approximation essentially involves accounting for the lift arising 
as a result of the geometric angle of attack at the tail associated with 
the pitching motion. 

The results also indicate that the simple time-lag effect, which is 
sometimes used to calculate the moment resulting from a constant vertical 
acceleration for aircraft at subsonic speeds, is not reliable at 

supersonic speeds. A more reliable method of calculating Cm. would be 

a 

to determine the upwash produced by a wing with a constant vertical 
acceleration. T hi s can be accomplished by the methods developed herein 
or by calculating the components of the upwash by the use of lifting 
lines. The contribution of the horizontal tail to Cm^ would then be 

found by the use of exact or approximate relations for the tail contri- 
bution to Cm* given herein. 


The results for two-dimensional wing-tail combinations seem to 
indicate that the damping of longitudinal oscillations due to C^ + 

is increased considerably if the tail surface lies between the Mach sheet 
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from the leading edge of the wing and the Mach sheet from the trailing 
edge of the wing. Calculated values of C m ^ + C ^ indicate that the 

damping contributed by this factor increases rapidly as the distance 
from the wing to the tail is increased. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , December 11, 1955* 
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APPENDIX A 


THE D 1 h'H'kkENTT ATTOW OF 02 


The function &£ is given by equation (22) or equation (2^). 
expression for 02 as given by equation (22) is 

0 2 (x,y,z) = — Z -~- f ■ J(x - J) 2 - B^Y 2 + zSjay-L + 


The 


2 nB 2 V^lii Y 2 + z 2 
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2 * 'Mttng (Y 2 + z 2 )R 


dxi dy 1 


By applying the rule for the differentiation of a definite integral with 
respect to a parameter to the first integral on the right-hand side of 
502 5©2 

equation (22), and can be written as 
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Equations (Al) and (A2) are valid when the limits of integration, bp 
and h 2 , are either the Mach cone from the point x,y,z or the plan-form 

tips. 


Equations (Al) and (A2) "become indeterminate at the apex of the 

hyperbola as z approaches zero. The value of ^ can 

directly from 0 2 (x,y,O) since this value is known. The value of 
de 2 (x,y,o) 


Sz 


is given by (from eq. (l7)> ref. 15) 


de 2 (x,y,o) 

Sz 



^ 2 A0 2 
dxp ctyi 
XT 


dx x dy-L - 


_1_ 

2rt 


d A9 2 

dx^ 

XY 


- dxp 


- £ A 
2 Sx 




Separating the wing and wake integrals yields 


d9 2 (x,y,0) 

bz 




P d £Q 2 


1_ 

2jt 


u) dxp 


d9. 


dxi - ® A 
XY 1 2 dx 


The second integral on the right-hand side of the preceding expression can 
be integrated with respect to xj_. The result of performing this integra- 
tion can be expressed as 
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d8 2 (x,y,0) 

dz 



Wing 


6x1 5 y x 
XY 


dx x dy x - 


1_ 

2 *, 


.2 

d A9 2 

\J{x - J) 2 - B 2 Y 2 

6x1 Sy-L 

Y 


- B cos 


-1 BY 


x - J 


<3yi 


Note that 


fi S A9p 
— R dxi 

lO ^ X 1 B A ^ 0 2 

2^ " 


XY 


- — A 


2 Sx 


(A3) 


5 A9 2 


may 


"be written as a function of y^_ since x^_ = J(yp). 


% X \ 6yi 

The expression for 9 2 as given by equation {2h) is 

0 2 (x,y,z) = f %(A,y,z)dA + -~ [^ <V(A,y,z)dA 
B V J M S B 2 V J Mg 


The lower limit of the preceding integrals is a function of y and z. 

In the differentiation of the preceding integrals with respect to either 
y or z, this variable limit must be considered. 

The trailing edge of the wing may be divided into three types. The 
first of these is the completely supersonic trailing edge. In this type 
of trailing edge, the component of free-stream velocity normal to the 
trailing edge is always supersonic. (See fig. 73(a).) The second type 
of trailing edge is the completely subsonic edge. Here, the component 
of free-stream velocity normal to the trailing edge is subsonic. (See 
fig. 75(b).) The third type of trailing edge is the mixed supersonic and 
subsonic trailing edge. As its name implies, this type of trailing edge • 
has both supersonic and subsonic portions. (See figs. 73(c) and 73(d).) 

For supersonic trailing edges the wake has no effect upstream of the 
trailing edgej thus, fip(x,y,z) is zero and equation (2k) reduces to 


0 2 (x,y,z) = -i- I n„.(A,y,z)dA {kh) 

B^V J Mg 
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In order to differentiate the integral of equation (A4 ) with respect to 
y and z it is convenient to express the lower limit Mg of the integral 

mathematically. For all hut the simplest trailing edges, the math ema tical 
difficulties are too great to make this expression practical. A general 
procedure is set up, however, and a few simple cases are presented. The 
equation of the trailing edge has been represented by 


X 1 = J(yi) 


and the Mach cone from each point on the trail ing edge is given by 


(x - J) 2 - B% 2 - B 2 z 2 = 0 


(A5) 


The Mach sheet from the trailing edge is the envelope of the Mach cones 
from all the points along the trailing edge. The equation of the Mach 
sheet from the trailing edge can be found by eliminating y^ from the 
two equations 


(x - J) 2 - B 2 (y - yq) 2 - B 2 z 2 = 0 (a6) 


a 

dyq 


(x - J) 2 



(A7) 


Such relationships follow from the mathematical procedure ' for finding the 
envelope of surfaces. (See ref. 23, p. 55-) 

When the trailing edge is swept at a constant angle and the t railing- 
edge end points do not affect the area of the Mach surface being con- 
sidered, the Mach surface is made up of two plane surfaces . For a 
coordinate system whose origin is located on the trailing edge, the 
equation of the trailing edge is given by 



Thus, 


J(yi) 


yi 

“t 
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and 

dJ(yi) = _1_ 

8yi “4: 

Equations (A 6) and (A 7 ) become 

( X " “ B2 ^ y " y l^ “ B2z2 = ° 



- B 2 (y - yq) = 


0 


The el imin ation of y-]_ from the two preceding equations yields 

(m-fcX - y) 2 - z 2 (B 2 mt 2 - 1) = 0 


Thus , the equation for the Mach surface is 


mt x - y - 



(A8) 


when z is positive and 


mtx - y + z 



0 


(A9) 


when z is negative. When the trailing edge is perpendicular to the 
stream direction, mq = «> and equation (a 8) becomes 


x - Bz = 0 


(A10) 


When the trailing edge is made up of a broken line composed of tWo 
straight-line segments, the Mach surface from the trailing edge is made 
up of parts of the after cones from the trailing-edge end points, parts 
of the after cone from the point connecting the two line segments, and 
the two Mach sheets from the two straight-line segments of the trailing 
edge. (See figs. 7^ and 75-) 
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In the preceding paragraphs a method has been presented for the 
determination of the equation of the Mach surface from the trailing 
edge. Let the equation of the Mach surface be represented by 


x = G(y,z) 

Equation (a 4) then can be written as 


0 2 (x,y,z) = - 5 - / %(\,y,z)dA 

B 2 V J G(y,z) 


(All) 


Differentiating equation (All) with respect to y and z yields 


de 2 (x,y,z) 1 [* ^r(\Y,z) _ 1 


dy -b 2 vJo(j,z) 


a* - sv(Ms)Mn5) 


B 2 V 


dy 


(A12) 


58 2 (x,y,z) = _p_ r* ^ 

Sz B 2 V^'G(y,z) ^z B^V 


» - -s- 


(A13) 


where fi^Mg ) is the value of ftw(^>y> z ) at the trailing-edge Mach sur- 
face obtained by approaching the Mach surface from the positive x-direction 
along the line (yp = y, zp = z). The expression iV(M 8 ) can be evaluated 

by calculating ft„. on the line (yp = y, zp = zj a small distance e 


downstream of the trailing-edge Mach 'sheet and then taking the limi t as 
e approaches zero. For points behind the Mach surface from any straight- 
line segment of a trailing edge, the quantity ft^Mg) is given by 




for positive z, provided the Mach surface at the point 


being considered is not part of the Mach cone from an end point of the 
line segment. 


For wings with subsonic trailing edges, ftp and ft„. at the trailing- 

edge Mach surface are zero. Equation ( 2k ) then can be differentiated to 
yield 
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For wings with m±xed supersonic and subsonic trailing edges, both 
Up and may be discontinuous through the Mach surface from the 


supersonic portions of the trailing edge. 


^0p 

Expressions for --=■ and 

oy 



for wings with mixed, trailing edges can be found in a manner similar to 
those used in treating wings with other types of trailing edges. The 
resulting expressions appear to be of little value since they require a 
knowledge of the induced potential on the plan form and for this reason 
they are not presented here. 
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APPENDIX B 

S0 Q _____ 

THE UPWASH — — ALONG THE WAKE ckntkk LINE OF A 

Bz 

PITCHING TRIANGULAR WING 


The upwash along the center line of the wake of a pitching triangu- 
lar wing with subsonic leading edges can "be found "by using a number of 
methods. The method used here is the potential doublet method presented 
in reference 11.' This method was chosen Because certain integrals which 
anise from using this method have already "been evaluated in reference 11. 

The potential difference across the surface of a triangular wing 
pitching about its apex is (from ref. 2*0 

AB q( x l-> y l) = 2riKqXi\]nPxi 2 ~ y± 2 (Bl) 


where the coordinate axes are located at the wing apex. (See fig. j6.) 

The potential at any point in the vertical plane of symmetry can be 
expressed as (ref. ll) 


0 a (x,O,z) = -B^zKq 


x l\j 


'bPxjZ - y x 2 


Wing 


[ 7 * - xp) 2 - B 2 (yp 2 + z2)j 


3/2 


dx-i d yi - 
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Equation (B2) can also "be expressed as 



The upwash along the center line of the wake is 

0 q 

w(x,0,0) = lim — 
z — »0 dz 
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The first integral of equation (B^) can be evaluated to yield 


1 1m 

z— >0 


TP — ii* U! 1 ) r * 1 ~ 71 ^ ^ ^ f ,v x . p *M - E (*a) 

Ju Wing ^ . ^2 . £^2 + i2 jj 3 /2 |_° (Bm + ij.) 2 Jl ^(toka + 

for c ^ x "k (Bm + l)c, and 


lla 
I“>0 


for (Bm + l)c = x. The last two integrals in equation (B*0 arise in 
evaluating the potential in the vertical plane of symmetry from a tri- 
angular wing. These two integrals were evaluated in reference 11. From 
reference 11, the values of these expressions are found to be given by 


-ltn B 2 *K a — 
x->0 ' d* 


X'fzPxjZ - 


yi* 


“•« \(* - *i) 2 - ^ 


3/2 


dy r - -2SqX 


P 1 *frl) - E ( k l) ^ . p ffe) - EfeQ : 

JO ^ + Bo Jl + Embj>) 


(B7) 


for c = x = (Bm + l)c, 

f, 


-lim 
z— >0 


zym^x-, 2 - yj 2 


B^A rr 

^Jjw^ j( x - x x ) 2 - B 2 ^ + Z 2 )J 




^ K ( kl ) ~ E ( k l) 


+ Bm 


<Ji x (b3) 


for (Bm + l) c ^ x. 


-lim Ak. A 

N.rt * OZ 


z — >0 


7= Z ^ 2 ~ yi2 -p Ti - - ^S-^3) - i 1 - k 3 2 ) K ( k 3]l 

( x - x l) 2 - B2( yi 2 + Z 2)j L - 1 


Wake 


(B9) 


for c i x k (Bm + l)c, and 


p 3 
-lim B^clC — 

z— ^0 dz 


z]jm 2 c 2 - y^ 2 


Wake 


{ x _ x j2 _ B 2( yi 2 + z 2) 


— dx-L dy-L = -2^01 (k 3 ) 


(BIO) 


for (Bm + l)c ^ x. 
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The upwash along the wake center line is (from eqs . (b 4) to (BIO)) 




for (Bm + l)c = x. 
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APPENDIX C 

THE UPWASH COMPONENT ALONG THE CENTER LINE OF 

dz 

THE WAKE OF A TRIANGULAR WING 


The upwash component 


^02 

dz 


can "be determined "by the methods pre- 


viously presented. For this case, however, the upwash can be determined 
relatively easily by the use of potential doublets and therefore this 
method is used. 


The difference in the x-derivative of the potential across the wake 
is given by 


A 



2a - y-^ 2 

B 2 E ' (Bm) 


(Cl) 


The potential difference across the wake is 




dx = 


2d(xi - c) ^n^e 2 - y-j 2 


B 2 E ' (Bm) 


(C2) 


The upwash component 
given by 


c)9p 

— =■ along the center line of the wake is 
oz 




66 


NACA TN 3072 


If the indicated operations in equation (C3) are performed, the following 
equations result: 


Se 2 

6z 


2oanc 
mBE ' (Bm) 




for c ^ x ^ (Bm + l)c, and 


60 2 
6z 


2a(x - c) 
mB 2 E ' (Bm) 


28 W 


- (1 - \ 2 ) k (\) 


for (Bm + l)c 6 x. 


( 0 +) 


(C5) 
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APPENDIX D 

EXPRESSIONS FOR THE MOMENT ON CERTAIN HORIZONTAL-TAIL SURFACES 


If the typ e of arguments utilized in references 22 and 25 are fol- 
lowed, expressions may he found for the lift and moment (about the trailing 
edge) on surfaces which have supersonic leading edges and trailing edges 
which are perpendicular to the free-stream direction (see fig. 77(a)) for 
the case of an angle-of -attack distribution which varies linearly with 
time. 


Consider a two-dimensional airfoil (as shown in fig. 78(a)) which 
has an angle of attack given by 


cr = -cr(t - to)h(x) 


(Dl) 


By symmetry it is clear that the lift and moment coefficients based on the 
area with the variable angle of attack for the wings shown in figures 78(a), 
78(b), and 78(c) are the same. The wing shown in figure 78(d) can be 
obtained as a simple combination of the wings shown in figures 78(a), 

78(b), and 78(c). The lift and moment coefficients based on the area 
with the variable angle of attack therefore are the same as the lift and 
moment coefficients for the two-dimensional wing. 

The potential function for the flow over the upper surface of the 
two-dimensional wing illustrated in figure 78(a) is (from eq. (13)) 


0 


r x_Bz 

b 3 ^0 


|h(|)d| + 


(t - t D ) 


M^x 

B 2 V 



fc(s)as 


(D2) 


where x is the distance downstream from the leading edge of the airfoil. 
The pressure-difference coefficient is 



[cr(t - t D )h(x) 


4-f 

B v v -'0 


h(|)d| 


(D3) 
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The lift and moment coefficients are, for t = to> 


C L - 


Lift /Unit length ter 


- V 2 ^ 


B^cV 


p C pX 

/ dx / h(|)d| (Ik) 

Jo J 0 




Moment /Unit length _ Ucr 


§V2S2 


• pC pX 

— - / (c - x)dx / h(g)dg (D5) 

’f v J 0 ^0 


3 _ 
Be V 


where the moment is taken about the trailing edge of the airfoil. 

The lift produced by a deflected strip such as is illustrated in 
figure 77(b) is the same as that produced by the strip illustrated in 
figure 78 (d-) and hence the same as a strip of the same width on the two- 
dimensional airfoil illustrated in figure J8 (a). The total lift and 
moment coefficients for the type of airfoil illustrated in figure 77(a) 
therefore are 


c L = - 


kb 


C m = 


ka 


-If 

[f 

b(6)ai c 

y Airfoil surface 


— 

JJ 

If 

h(g)dg 

Airfoil surface 

J L.E. 



ds 


ds 


(D6) 


(D7) 


where A is the distance upstream to the element of area ds. (See 
fig. 77(b).) 

In the case where the airfoil illustrated, in figure 77 ( a ) is a tail 
surface and the wing-tail combination has a constant vertical accelera- 
tion, the time -dependent angle-of -attack distribution on the tail surface 
is given by 


a 


a 1 + 


8 fia' 
9z aV, 


t 
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In this case equations (l>6) and (D 7 ) become 


c L 


Lift 
| V 2 S t 


— ff 

B^S-j-V Tail surface 





ds 


(D8) 


^ = 


Moment 

§ v2s #t 


4a 


B 5 S t c t V 


Tail surface 




+ 


-—Vi 

dz aV/ 


ds 


(D9) 
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TABLE I.- F0R4DLAS FOR TRANSFORMATION OF THE STABILITY DERIVATIVES C^, 
C]^, C^, C]- , AND FROM BODY TO STABILITY AXES 


Body axes; 
origin^ a t 
x - 0, y = 0, 
z » 0 

at 

Body axes; origin 
x = d, y « 0, z = 0 


Stability axes; origin 
at x *» d, y = 0_, z - 0 

Stability 

derivative 

Stability 

derivative 2 

Shift in origin 
from (0,0j0) to (d,0,0) 

Stability 

derivative? 

Origin at (d.,0,0); rotation 
through angle a 


V 

^ + 1“ °Iu 
Cv 


V 

C^* (approx.) 

% 

<v 



V 

Cl^* (approx.) 

% 

V 


°Iu 

S' 

V 

°Id 


' °id 



Cj^* (approx.) 

Sx 




<V 

<*■£* 


^The origin refers to a system of Cartesian coordinates used in the analysis. (See fig. 25 (a)*) 

2 These coefficients are for a system of Body axes located at (d,0,0) in the system used in 
the analysis. (See fig. 25(a).) 

^These coefficients are for a system of stability axes located at (d,0,0) in the 
system used in the analysis. (See fig. 25(b).) 
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Figure 8.- Exact (linearized) and approximate (eq.. (3Tb)) values of 

upwash at t = t Q along the center line of the wake behind a triangular 
wing with constant vertical acceleration for various Mach numbers. A = 2. 
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Figure 9-- Location of coordinate axes used in treating a rectangular 

vine tin. 
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(d) AB = 3. 


Figure 11.- Continued. 
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(e) AB = k. 
Figure 11.- Continued. 
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b/Z 


(f) AB = 6. 


Figure 11.- Continued. 
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(g) AB = 8. 
Figure 11.- Concluded. 
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(a) x = 2c; M = 1.20. 

Figure 16,- The contribution of each component of upwash in the plane of 
the wing at t = t Q for a eemi-infinite rectangular wing with a con- 
stant vertical acceleration. 
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(c) x=«2c;M = 1.4l4. 

Figure 1 6 .- Continued. 
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(d) x = 4c; M = 1.414 

Figure 16.- Continued. 
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(a) M = 1.414. 


Figure 17 .- The upwash distribution at t = t Q in the plane of the wing 

be hin d a rect an g u l ar wing with a constant vertical acceleration with 
A = 1 for various Mach numbers. 
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(b) M = I.667. 


Figure 17-- Continued. 
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(c) M = 2.24. 


Figure 17.- Concluded. 



n8 


MCA TN 3072 



Figure 18 .- The upwash distrihution at t = t Q in the plane of the wing 
hehind a rectangular wing with a constant vertical acceleration with 
A = 2 for various Mach numbers . 
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Ob) M = l.klk. 


Figure 18.- Continued. 









(a) M = 1.20. 


Figure 19- - The upwash distribution at t = t Q in the plane of the wing 

behind a rectangular wing with a constant vertical acceleration with 
A = 3 for various Mach numbers. 
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(a) M = 1.20. 

Figure 20.- The upwash distribution at t = t Q in the plane of the wing 
behind a rectangular wing with a constant vertical acceleration with 
A = 4 for various Mach numbers. 
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Figure 21.- The upwash distribution at t = t Q in the plane of the wing 
behind a rectangular wing with a constant vertical acceleration with 
A = 6 for various Mach numbers. 




I 


l 

1 


1 

1 



I 



ft) M - 1*^* 

2i,- Continued. 


& 

l - 1 


Figure 


sa 





(a) M = 1.20. 


Figure 22.- Tbe upwash distribution at t = t Q in the plane of the wing 
behind a rectangular wing with a constant vertical acceleration with 
A = 8 for various Mach numbers. 
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(a) M = 1.4-14- . 

Figure 23.- Exact (linearized) and approximate (eq. (37^)) values of 
upwash at t = t Q in the plane of the wing behind a rectangular 
wing for two Mach numbers. A = 1. 
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(a) M = 1.414. 

Figure 24.- Exact (linearized) and approx ima te (eq. 
upwash at t ■= t G in the plane of the wing hehi 
for two Mach numbers. A t 3. 
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(b) Stability axes (principal body axes dashed for comparison) . 
Figure 25.- Systems of axes and associated data. 
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(a) Tail lies downstream of Mach sheet from trailing edge of wing. 




(t>) Mach sheet from wing trailing edge intersects tail. 



(c) Tail surface lies between Mach sheet from leading edge of wing and 
Mach sheet from trailing edge of wing. 


Figure 2 6 .- Types of flow treated for the two-dimensional 
wing-tail combination. 




,ed) and approximate (eq.. (55)) 
wing-tall combination, at two 












c w 

Figure JO.- Comparison of exact (linearized) and approximate (eq. (53)) 
values of for a triangular wing-tail combination at two 

Mach numbers. 
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Wing A = / 
Tail A - 4 



Wing A = 3 
Tail A - 4 
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Wing A - 2 
Tail A - 4 



Wing A - 4 
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Figure '51.- Several re et ang ul ar -wing — triangular-tail combinations and 

associated data. 



























(a) Wing aspect ratj 


Figure 3^.- Comparison of exact 
values of C™. for two rect 



$ Approx., M= 1.41 

a Approx., M=3.!6 


M = 3.16 
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Exact 

Approx v M * 1,41 
Approx M*2.24 








Figure 35*- Comparison of exact (linearized) and approximate (eq.. ( 65 )) 
values of C^ for a tvo-dimensional wing-tall combination for two 

Ctj- 

Maci. numbers. X = c+ = 
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Figure 56.- Comparison of exact (linearized) and approximate (eq.. (66)) 
values of for a tvo-dimensional wing-tail combination for two 

Macb numbers . x = 




-24 




Figure 58.- Variation of 









Figure 59.- Comparison of exact (linearized) and approximate (eq.-. ( 65 )) 
values of Cr for a triangular wing-tail combination at two Mach 
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Figure ^0.- Comparison of exact (linearized) and approximate (eq.. (66)) 
values of Cm^ for a triangular wing -tail combination at two Mach 


numbers . 
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Figure 43.- Comparison of exact (linearized) and approximate (eq_. ( 65 )) 
values of for two center-of -gravity locations for a triangular 

wing-tail combination. M = 1.1)4. 
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(a) Wing aspect ratio 1; — = i. 

Sr 2 

Figure 47-- Comparison of exact (l i nearized) and approximate (eg.. ( 65 )) 
values of Cp, for two rectangular-wing-- triangular-tail combinations. 
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(a) Wing aspect ratio 1; — = — . 

°v 2 

Figure ^9*- Variation of with l/cy. for three center-of -gravity 

locations for two rectangular-wing-*. triangular -tail combinations. 
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(t) Wing aspect ratio 4; — = 1. 
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Figure 50.~ Concluded. 
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(a) Wing aspect ratio 1 : — a i. 
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Figure 51 .- Comparison of exact (linearized) and appro xima te (eq.. (65)) 
values of Cl^ for two center-of -gravity locations for two rectangular - 

wing--trlangular-tail combinations. M o l.lt-11).. 
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(a) Wing aspect ratio 1; — t= — . 

°v 2 

Figure 52.- Comparison of exact (linearized) and approximate (eq. (66)) 
values of for two center-of-gravity locations for two rectangular- 

wing — triangular -tail combinations. M = 1.4l4. 







Figure 53 •- Comparison of exact (linearized) and approximate (eq. (8l)) 
values of Cp^ for a two-dimensional wing-tail combination at two 

, - _ — Cxr 

Macn manners . x = c-fc = 
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Figure 55.- Variation of with l/c^ for three Mach numbers for a 

triangular wing -tail combination. 
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Figure 57 •- Comparison of exact (linearized) and approximate (eq.. (8l)) 

values of Cl^ for a triangular wing-tail combination at two Mach H 

numbers . 







Figure 59.- Variation of Cm^ with l/c^ for three center-of -gravity 
locations for a triangular wing-tail combination. M = 1.44. 
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Figure bO.- Comparison between exact (linearized) anrl approximate 
(eq. ( 82 )) values of for two center-of-gravity locations for 

a triangular wing-tail combination. M = 1.44. 
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(c) Wing aspect ratio 2j = 1. 
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Figure 6l.~ Continued. 
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(a) Wing aspect ratio lj — « — . 

% 2 

Figure 6k. - Comparison between exact (linearized) and approximate (eg. (82)) 
values of C^ for two rectangular -wing -^triangular-tail combinations. 
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Figure 65 .- Variation of with l/c^. for three center-of -gravity 

locations for two rectangular -wing — triangular-tail combinations. § 
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(a) Wing aspect ratio lj ^ a i. 
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Figure 66.- Comparison between exact (linearized) and approximate (eq.. (82)) 
values of for two center-of-gravity locations for two rectangular- 
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Figure 68.- Variation of 




ifcy for a two-dimensional 


ving-tail combination for two Mach, numbers. 
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Figure 70.- Variation of Mach number M 
tion x/c^. for which Cm^ + Cm^ = 0 
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with center -of -gravity loca- 
tor five two-dimensional- wing- 


tail combinations 
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(a) Plan form with a supersonic trailing edge. 



(b) Plan form with a subsonic trailing edge. 



(c) Plan form with mixed supersonic and subsonic tra i l in g edge. 



(d) Plan form with mixed supersonic and subsonic trailing edge. 
Figure 73.- Plan forms with various types of trailing edges. 
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(a) Arbitrary angle -of -attack distribution over complete plan form. 



er-o-h(X)(t-L) 


(b) Arbitrary angle-of -attack distribution confined to a strip. 

Figure 77 •“ Airfoils with supersonic leading edges and trailing edges per 
pendicular to the free-stream direction. 
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(a) Arbitrary angle -of -attack distribution-. 



(b) Arbitrary angle -of -attack distribution confined to left side of airfoil. 



(c) Arbitrary angle-of -attack distribution confined to right side of airfoil. 



(d) Arbitrary angle-of -attack distribution confined to a strip. 

Figure 78 .- Infinite-aspect-ratio airfoils with arbitrary angle -of -attack 
distributions varying linearly with time. 
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